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ABSTRACT 

A  procedure  is  described  which  permits  the  computation 
of  the  radiation  characteristics  of  a  transducer  array  given  the 
electrical  inputs  to  the  transducer  elements.  The  effects  of 
transducer  element  interactions  (acoustic  coupling)  and  dome- 
transducer  interactions  are  included  in  the  computation. 

The  dome -transducer  configuration  chosen  for  this  study 
is  a  rectangular  array  of  circular  elements  mounted  in  a  spher¬ 
ical  baffle.  This  array  is  surrounded  by  a  concentric,  homog¬ 
eneous  dome.  Numerical  results  include  (a)  radiation  loadings 
on  the  transducer  elements,  (b)  response  (head  velocity)  of  each 
element  to  prescribed  electrical  inputs  and  (c)  farfield  beam 
patterns  for  the  dome/array.  The  results  indicate  that  element 
interactions  and  dome/array  interactions  can  affect  element  re¬ 
sponse  to  the  electrical  input  and  thereby  degrade  transmit  per¬ 
formance  of  the  transducer  array. 
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1.  INTRODUCTION 

The  effects  of  acoustic  coupling  among  transducer  ele¬ 
ments  on  the  performance  of  large  sonar  transducer  arrays  oper¬ 
ating  at  low  frequencies  have  been  the  subject  of  extensive 
research.  In  independent  studies,  the  effect  of  a  dome  on  near¬ 
field  and  farfield  radiation  characteristics  of  such  arrays  has 
been  investigated  assuming  negligible  coupling  effects  among  the 
elements.  It  has  been  found  that  the  presence  of  a  dome  can 
significantly  alter  the  radiation  loading  on  an  active  array 
relative  to  the  loading  in  the  absence  of  the  dome.  The  purpose 
of  this  memorandum  is  to  illustrate  a  technique  for  determining 
the  effects  of  acoustic  coupling,  both  among  transducer  elements 
and  between  the  dome  and  array,  on  array  performance. 

When  a  single  transducer  element  is  radiating  into  a 
free  field,  the  radiation  impedance  seen  by  the  element  is  essen¬ 
tially  pc,  the  characteristic  impedance  of  the  surrounding  medium. 
If  a  specified  array  of  elements  is  considered,  the  radiation 
loading  on  any  one  element  depends  on  several  factors,  including 
element  location  in  the  array,  array  configuration  (cylindrical, 
spherical,  etc.),  beam  steering/tilt  angle,  and  frequency.  When 
the  array  is  surrounded  by  a  dome,  the  element  radiation  loadings 
are  modified  further.  The  radiation  loading  on  any  single  ele¬ 
ment  can  determine  the  response  (head  velocity)  of  that  element, 
assuming  a  known  driving  voltage  or  current.  Since  the  element 
loadings  can  vary  considerably  in  the  array,  the  velocity  dis¬ 
tribution  of  the  array  and,  consequently,  the  radiation  charac¬ 
teristics  of  the  dome/array  can  be  affected  compared  to  the  ideal 
case  for  which  an  element  response  is  uniquely  determined  by  that 
element's  applied  electrical  signal. 

To  analyze  element  response  in  the  dome/array  environ¬ 
ment  and  subsequently  compute  dome/array  radiation  characteris¬ 
tics,  one  must  have  (1)  a  dome/array  model  which  can  treat 
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radiation  due  to  an  arbitrary  number  of  elements  having  arbitrary 
head  velocities,  and  (2)  a  transducer  element  model  which  yields 
response  as  a  function  of  radiation  loading  for  prescribed  input 
voltage  or  current.  The  procedure  is  demonstrated  in  this  memo¬ 
randum  through  utilization  of  a  spherical  array  surrounded  by  a 
concentric  shell  and  an  assumed  transducer  element  model. 
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2 .  BACKGROUND 


Transducer  element  modeling  has  progressed  in  stages 
from  crude  lumped  parameter  models  to  more  sophisticated  ones 
in  which  components  of  the  element  arc  treated  as  continuous 
media*.  As  one  goal,  these  models  seek  to  relate  realistically 
the  response  of  an  element  to  a  prescribed  electrical  signal  as 
a  function  of  radiation  loading.  Often  these  models  are  repre¬ 
sented  by  equivalent  circuits.  In  this  form,  the  input  variables, 
voltage  and  current,  can  be  related  to  the  output  variables,  head 
velocity  and  radiation  force,  by  a  pair  of  linear  independent 
equations.  The  coefficients  in  these  equations  are  obtained 
through  mathematical  and  experimental  analysis  of  the  element. 

The  element  modeling  has  been  accompanied  by  much 
analytical  and  experimental  effort  to  determine  the  effect  of 
the  array  environment  (radiation  loading)  on  element  response**. 

If  we  consider  an  array  of  elements  and  assume  the  ith  element 
has  a  nead  velocity  v.  and  as  a  result  produces  a  sound  pressure 
field  p^,  then  a  force  is  created  on  each  of  the  neighboring 
elements.  In  particular,  the  force  F..  on  the  jth  element  is 


p .  dS . 
J 


(1) 


•*>  S .  is  the  area  of  the  jth  element.  If  F..  is  normalized 
J  J  iJ 

then  we  define  the  interaction  coefficient  Z. .  between  the 

ith  and  jth  pistons  as 


when 
by  v± 


*D .  G.  Carson,  G.  E.  Martin,  et.  al. ,  of  the  Transducer  Division, 
NUWC  have  been  principal  contributors  in  this  area. 

**Among  others,  the  NUWC  group  and  D.  T.  Porter  of  USN/USL  have 
been  active  in  this  area. 
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Z.  . 
1J 


(2) 


with  Z.  =  Z. . . 
i.;i  .]  i 

The  total  force  F.  on  the  ith  element  due  to  N  active 

J  J 

elements  in  the  array  is* 

N  N 


F. 

J 


i=l 


i=l 


with  Z..  v.  the  radiation  force  on  the  ith  element  due  to  its 
JJ  J  rad 

own  motion.  The  total  radiation  impedance  Z.  seen  by  the  jth 

element  is  given  as 


zfad 

J 


(4) 


If  the  internal  impedance  of  the  element  is  Zqc  (the  same  for  all 
elements  in  the  array) ,  then  the  head  velocity  v^  can  be  related 
(see  Ref.  1)  to  the  current  1^  (or  voltage)  applied  to  the  ith 
element  by 

I.  =Y(Z^d  +  Zoc)v.  (5) 


where  y  is  a  complex  constant  dependent  on  the  element  component 
materials,  dimensions,  static  biases,  etc.,  and  the  operating 
frequency.  Note  from  Eq.  (5)  that,  in  general,  the  head  velocity 

of  the  ith  element  depends  on  the  head  velocities  of  all  the  ele¬ 
ments,  i.e.,  the  currents  and  head  velocities  are  related  by  a 
set  of  simultaneous  equations.  In  matrix  notation  Eq.  (5)  is 
represented  as 

'^Interaction  effects  due  to  nonactive  elements  are  not  considered 
here,  although  they  can  be  included  in  the  analysis. 
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{  1  il  =  [ct^]  £V-jd  >  (6) 

with  a*.  =  YZ..,  i^j  and  a-.  =  y (Z . .  +  Z  ).  A  similar  relation 

J  JL  I  X.  JL  XX  U  X- 

exists  for  voltage  controlled  elements. 

The  above  relation  permits  the  computation  of  a  set  of 

head  velocities  for  a  prescribed  set  of  currents,  assuming  the 

Z..  can  be  computed.  For  the  computation  of  the  Z..,  a  mathemati- 
— *  J 

cal  model  is  required  which  simulates  the  dome/array  of  interest. 

The  array  model  must  allow  the  computation  of  the  sound  pressure 
produced  at  any  point  on  the  array  by  an  element  having  arbitrary 
location  in  the  array  and  arbitrary  head  velocity.  This  model 
is  exercised  between  all  element  pairs  to  obtain  the  set  of  Z... 
(Symmetry  considerations  often  allow  a  reduction  in  the  required 
number  of  computations.)  A  set  of  head  velocities  is  then  deter¬ 
mined  from  Eq.  (6)  for  known  {1^},  Zqc,  and  y.  This  set  of  head 
velocities  is  utilized  in  the  dome/array  model  to  compute  near¬ 
field  and  farfield  sound  pressure  levels,  beam  patterns,  etc. 
Descriptions  of  several  dome/array  models  which  can  be  used  in 
such  a  series  of  computations  are  given  in  Ref.  2. 
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3.  TECHNICAL  DESCRIPTION 

The  dome/array  model  chosen  for  this  study  consists  of 
a  rectangular  array  of  close-packed,  circular  elements  on  a 
spherical  baffle.  The  baffle  is  surrounded  by  a  concentric  shell 
as  shown  in  Fig.  1.  Radiation  characteristics  for  this  model, 
assuming  element  velocities  independent  of  radiation  loading, 
have  been  studied  in  Ref.  3. 

The  analysis  is  initiated  by  considering  the  sound 
pressure  field  produced  by  a  single  circular  piston  or  element 
radiating  with  harmonic  time  dependence.  As  shown  in  Ref.  3, 
the  sound  pressure  field  p^  between  the  array  and  dome  is  given 
by 


pi 


CO 


m=o 


[A  j  (kr)  + 
mJm '  ' 


B  n  (kr)] 
m  mv  7 J 


P  (cos  \[i )  e 
m 


-iwt 


(7) 


and  the  sound  pressure  exterior  to  the  dome,  p£,  is  given  by 

00 

p2  =  1  Cmhti/kr)  Pit/C0S  e~mt  ‘  (8) 

m=o 


In  the  above  expressions,  jm(kr)  ,  nm(kr) ,  and  hm(kr)  are  the  mth 
order  spherical  Bessel,  Neumann,  and  Hankel  functions  respectively, 
k  is  the  wave  number,  P  (cos  ty )  is  the  mth  order  Legendre  poly¬ 
nomial,  is  the  polar  angle  measured  from  an  axis  through  the 
center  of  the  element  (the  problem  is  axisymmetric) ,  and  u)  is 
the  angular  frequency.  A  ,  B  ,  and  are  constants  determined 
by  the  boundary  conditions.  The  boundary  conditions  are  defined 
by  specifying  the  normal  particle  velocity  of  the  fluid  to  be 
zero  everywhere  on  the  sphere,  except  for  the  piston,  which  has 
a  uniform  velocity  Voe*"la't,  and  requiring  the  normal  particle 
velocity  at  fluid-shell  interfaces  to  be  continuous. 
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Three  basic  steps  are  involved  in  obtaining  radiation 
characteristics  of  the  dome/array  from  the  known  electrical  inputs 
to  the  transducer  elements.  These  are: 

1.  The  loading  produced  by  any  element  on  the  remaining 
elements  in  the  array  is  obtained  by  integrating 

Eq.  (7),  evaluated  on  the  array  surface,  over  each 
element  in  the  array  for  all  pairs  of  elements  as 
in  Eq.  (1) .  These  loadings  are  then  converted  to 
as  in  Eq .  (2)  . 

2.  The  set  of  head  velocities  is  obtained  by  solving 
the  set  of  equations  given  by  Eq.  (5)  for  known 
input  currents  (or  voltages) ,  element  characteris¬ 
tics,  and  the  Z..  obtained  under  (1)  above. 

’  ij 

3.  Equation  (7)  is  evaluated  for  each  element  using 
the  head  velocity  computed  in  (2)  above,  and  the 
resulting  element  pressure  fields  are  superimposed 
to  obtain  the  sound  pressure  field  for  the  array 
between  the  array  and  dome.  The  sound  pressure 
field  exterior  to  the  dome  is  obtained  as  a  super¬ 
position  of  the  individual  element  pressures  given 
by  Eq.  (8)  . 

The  first  step,  determining  the  Z„ ,  is  straightforward. 
As  shown  in  the  appendix,  the  sound  pressure  field  produced  by  a 
single  piston  has  the  form 

P1  =  Vo  _  rn!"Pn-l^cos  a)  “  Pn+l^cos  a^]  Pn^cos  ^  ^ 

n 

when  evaluated  on  the  array  surface.  In  the  above  equation,  the 
rn  are  constants  depending  on  the  dome/array  dimensions,  dome 
material  properties,  and  frequency,  and  a  is  the  angular  half¬ 
width  of  the  element  (if  the  element  radius  is  a  and  the  array 

SLq  ® 

radius  is  a, then  a  =  arctan  — )  .  To  obtain  the  loading  on  the 

Cl 
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jth  element  in  the  array,  separated  from  the  ith  element  by  an 
angle  9 ,  we  integrate  Eq.  (9)  over  the  jth  element  and  obtain 


00 


F  .  .  =  V  n  a  o  c 
lj  o  „  n 

n=o 


[Pn_1(cos  a)  -  Pn+1(cos  a)] 

2_I 


Pn(cos  9).  (10) 


F..  will  be  normalized  by  ra  pc  so  that 
lj  oK 


2  “ 

Z.  .  =  K 

lj  2  __ 

J  a 

o  n=o 


n 


[Pn_1(cos  a)  -  Pn+1(cos  a)] 
2n+I 


Pn(cos  8)  .  (11) 


The  are  computed  for  all  pairs  of  active  elements  in  the 
array. 

The  next  step  in  the  series  of  computations  is  the 

determination  of  the  set  of  element  head  velocities.  For  these 

173(1 

computations,  Z  for  each  element  and  Zqc  must  be  determined. 

Zqc  is  obtained  from  the  transducer  element  model.  An  example 
of  the  type  of  element  considered  in  this  memorandum  is  shown 
in  Fig.  2.  Along  with  a  cross  section  view,  a  circuit  schematic 
is  shown.  The  components  are  assumed  to  be  one-dimensional  media. 
Wave  propagation  in  the  media  is  described  by  linear  elasticity 
theory.  The  components  of  the  element  must  satisfy  boundary 
conditions  at  their  juncture  with  other  components.  For  example, 
forces  and  velocities  must  be  continuous  at  the  boundaries  of 
each  component.  For  the  passive  components  (head  mass,  tail 
mass,  etc.),  displacements  satisfy  the  scalar  wave  equation,  but 
for  the  active  component  (ceramic) ,  a  mechanical  strain  component 
due  to  the  applied  electrical  signal  must  be  considered.  Solu¬ 
tions  to  the  governing  differential  equations  can  be  obtained  by 
usual  techniques.  Due  to  the  coupling  of  the  element  components 
at  the  boundaries  (e.g.  stress  rod  -  tail  mass  interface),  the 
equations  describing  component  motion  are  coupled.  The  derivation 
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of  the  element  model  and  an  accompanying  equivalent  circuit  are 
shown  in  detail  in  Ref.  1. 

For  this  study,  an  element  model  similar  to  the  one 
shown  in  Fig.  2  and  analyzed  in  Ref.  1  was  adopted.  This  ele¬ 
ment  has  a  Zqc  which  varies  with  frequency  as  shown  in  Fig.  3. 
(Zqc  is  plotted  as  a  function  of  ka,  where  k  is  the  wave  -.umber 
and  a  is  the  radius  of  the  sphere).  For  a  given  ka,  Zqc  is 
obtained  from  this  curve,  the  constant  Y  is  computed,  and  the 
Z^j  (also  a  function  of  ka)  are  computed.  These  quantities, 
together  with  the  prescribed  electrical  signals,  are  then  used 
to  compute  a  set  of  head  velocities  {v^}  from  Eq.  (6). 

The  final  step  in  the  computations  is  the  evaluation 
of  the  radiation  characteristics  of  the  array.  First,  the  pres¬ 
sure  fields  of  the  individual  elements  are  computed  from  Eq.  (7) 
(nearfield)  or  Eq.  (8)  (farfield) .  These  are  then  superimposed 
to  obtain  the  pressure  field  for  the  array.  A  simple  angle 
transformation  is  used  in  the  superposition,  as  discussed  in  the 
Appendix.  This  superposition  of  the  element  pressure  fields 
completes  the  series  of  computations. 
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4.  NUMERICAL  RESULTS 

Numerical  Results  are  presented  for  three  cases;  (1) 
at  the  "tuned"  frequency  fQ,  for  ka  =  36.0,  where  Zqc  is  a  maxi¬ 
mum,  (2)  at  a  frequency  f^,  below  fQ,  at  which  ka  =  31.9,  and 
(3)  for  "ideal"  head  velocities,  i.e.,  the  head  velocity  of  each 
element  is  affected  only  by  the  self-impedance  of  that  element. 
Results  are  shown  both  with  and  without  dome  for  each  of  the 
three  cases. 

The  transducer  used  in  the  computations  consists  of  an 
array  of  8  x  24  half-wavelength  (at  f  ,  ka  =  36.0),  close-packed 
circular  elements.  The  dome  has  kb  =  51.0  and  kh  =  0.1164  at  f 

o 

and  has  a  specific  gravity  of  7.86  (steel).  The  input  currents 
{ 1^}  are  of  equal  unit  amplitude  and  phase  delayed  to  a  plane 
for  each  frequency. 

Figures  4-7  show  contour  maps  of  the  magnitude  of  the 
total  radiation  loading  on  the  array.  The  data  in  each  curve 
are  normalized  within  that  data  set.  The  range  of  loading  from 
the  minimum  to  maximum  is  divided  into  equal  increments  with 
symbol  "l"  corresponding  to  the  lowest  value  and  "A"  corresponding 
to  the  highest  value.  The  effect  of  a  dome  on  radiation  loadings 
is  graphically  illustrated  by  a  comparison  of  Figs.  6  and  7.  For 
frequency  f^  and  no  dome,  the  radiation  loadings  are  nearly  con¬ 
stant  over  the  broad,  central  portion  of  the  array,  and  increase 
smoothly  to  maximum  values  at  the  corners.  However,  the  addition 
of  a  dome  (Fig.  7)  causes  a  large  buildup  of  radiation  loading 
in  the  central  portion  of  the  array.  Note  that  the  radiation 
loading  exceeds  Z  by  an  order  of  magnitude  for  the  case  shown 
in  Fig.  7. 

The  normalized  magnitudes  of  the  velocities  at  frequency 
fL>  with  a  dome,  for  two  rows  of  the  array  are  shown  in  Fig.  8. 
Comparing  Fig.  7  and  Fig.  8  shows  good  correlation  between  them, 
in  that  elements  with  less  radiation  loading  have  higher  velocity 
magnitudes . 
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The  farfield  beam  patterns  are  presented  in  Figs.  9-11 
with  patterns  for  both  ideal  and  computed  velocities  superimposed 
on  the  same  figure.  The  beam  patterns  are  normalized  and  do  not 
show  the  effects  of  acoustic  coupling  on  insertion  loss.  These 
patterns  show  that  including  element  interaction  results  here  in 
a  slightly  narrower  main  lobe  and  increased  side  lobe  levels. 

The  most  serious  difference  is  the  case  (Fig.  11)  with  frequency 
fT  ,  with  a  dome,  where  the  first  side  lobe  is  up  nearly  10  dB 
from  the  ideal  case,  and  the  main  lobe  is  about  2  narrower. 

Figure  12  shows  the  insertion  loss  (change  in  intensity  on  the 
beam  axis)  for  each  beam  pattern  computed  at  fQ  relative  to  the 
no-dome  case  with  "ideal"  velocities  at  fQ.  It  can  be  seen  that 
both  dome  interaction  and  element  interaction  influence  insertion 
loss.  The  relative  importance  of  dome-transducer  interactions 
and  element  interactions  in  determining  insertion  loss  should  not 
be  generalized  from  these  results  because  of  the  strong  dependence 
of  these  interactions  on  dome-transducer  geometry,  frequency,  and 
element  size  and  spacing. 
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5.  CONCLUSIONS 


This  memorandum  demonstrates  a  procedure  for  analyzing 
the  effects  of  transducer  element  interactions  and  dome/array 
interactions  on  element  response  and,  consequently,  radiation 
characteristics  of  the  array.  As  shown  in  the  computed  results, 
these  interactions  can  have  significant  effect  on  the  performance 
of  an  array.  Realistic  analysis  of  interaction  effects  is  par¬ 
ticularly  important  if  the  array  must  perform  over  a  frequency 
band.  In  this  case  the  internal  impedance  of  the  element, 
usually  a  maximum  over  a  narrow  range  of  frequency,  and  the 
radiation  impedance  can  be  of  similar  magnitude  at  the  upper  and 
lower  limits  of  the  operational  frequency  band.  As  a  result, 
the  radiation  characteristics  of  the  array  can  be  degraded. 
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APPENDIX 

SPHERICAL  TRANSDUCER-DOME  MODEL 

The  analysis  of  the  spherical  transducer  with  concentric 
spherical  dome,  outlined  here,  is  presented  in  a  previous  tech¬ 
nical  memorandum  [Ref.  3].  The  geometry  of  the  spherical  trans¬ 
ducer  and  dome  is  shown  in  Fig.  1. 

The  first  step  in  the  analysis  is  to  obtain  the  solution 
for  the  sound  pressure  radiated  by  a  single  circular  piston.  A 
piston  vibrating  with  harmonic  time  dependence  will  produce  a 
sound,  pressure  field  between  the  dome  and  transducer  which  can 

\  t  f- 

be  represented  by  the  expansion  (in  which  e~  U;  time  dependence 
is  understood) , 

co 

p.  =  V  ("a  j  (kr)  +  B  n  (kr)~|  P  (cos  4).  (A-l) 

rl  L  L  mJmv  '  mmx/Jmv  Y/  v  ' 

m=0 


This  form  is  standard  for  solutions  of  the  wave  equation  in  spheri¬ 
cal  coordinates  (r,  f)  if  there  is  symmetry  about  the  axis  from 
which  the  angle  ijj  is  measured.  The  P^  are  Legendre  polynomials 
of  argur --»nt  cos  and  the  jm  and  nm  are,  respectively,  the  spheri¬ 
cal  Bessel  and  Neumann  functions  of  argument  kr.  The  A  and  B 
are  constants  which  can  be  adjusted  to  make  the  solution  form 
fit  whatever  boundary  conditions  are  prescribed  on  the  transducer 
and  shell  surfaces.  The  boundary  conditions  arise  in  equating 
the  fluid  and  solid  velocity  components  normal  to  the  fluid- solid 
interfaces.  Velocity  components  tangential  to  the  interfaces  are 
not  necessarily  equal  since  the  fluid  is  assumed  inviscid. 

For  the  case  at  hand,  the  fluid  velocity  is  equal  to  the 
piston  velocity  on  the  piston  surface  and  is  zero  on  the  remain- 
ing  transducer  surface.  The  piston  velocity,  vQe  ,  will  be 
considered  to  be  purely  radial,  so  the  boundary  condition  at  the 
transducer  surface  r  =  a  is 


A-l 


v0,  0  S  ;■  <  'J 

<  .  (A-2) 

0,  a  <  \jr  £  n 

Here  the  angular  extent  of  the  piston  is  0  ^  \!»  <  a  .  the  radial 

fluid  velocity  component  being  zero  outside  of  this  range.  The 

primes  on  the  i  and  n  denote  differentiation  with  respect  to 
v  Jm  m  r 

argument . 

While  the  pressure  field  between  the  transducer  and 
dome  has  the  form  given  in  Eq.  (A-2),  the  sound  pressure  field, 

P£,  outside  the  dome  must  have  a  form  which  describes  only  out¬ 
going  radiation.  This  form  is 


p2  =  Cmhm^kr^  Pm^C0S  ^  ’  (A“3) 

m=0 


where  h  =i  +in  is  a  spherical  Hankel  function  and  the  C 
are  determined  by  the  boundary  condition  at  the  outer  dome  sur¬ 
face.  Since  the  dome  is  thin  relative  to  an  acoustic  wavelength, 
it  is  sufficient  to  apply  the  velocity  continuity  condition  at 
the  dome  midsurface,  r  =  b,  rather  than  at  the  actual  fluid- dome 
interfaces.  Thus, 


1  5pl  =  _1_  1^2  =  w 

imp  a  r  — b  imp  a  r  ^ 


(A-4) 


w  being  the  radial  velocity  distribution  of  the  dome.  From  the 
given  boundary  conditions  and  the  analysis  of  the  dome  motion 
given  in  the  earlier  technical  memorandum  rRef.  3]  the  unknown 
coefficients  Am,  Bm,  and  Cm  can  be  found.  Equations  (A-l)  and 
(A-3)  can  now  be  used  to  find  nearfield  and  farfield  pressures 
respectively . 

The  pressure  field,  p^,  generated  by  a  single  piston 
and  evaluated  at  the  transducer  surface,  is  found  to  be 
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?!  "  Vo 


rn[pn-l(cos  a)  -  Pn+1(cos  a)]  Pn(cos  i })  ,  (A- 5) 


where 


hn(ka)  +  -^-(kb)Zhn(kb)[  jn(ka)n^(kb) 

Z  ~ 

h’  (ka)  +  t— —  (kb)  ^h!  (kb)[  j  1  (ka)n'  (kb) 
n v  7  ipcx  7  n  /UJnv  7  n  7 


j^(kb)nn(ka)] 

jn(kb)nn(ka)] 


and  is  the  modal  impedance  of  the  shell  (see  Ref.  3) .  This 
pressure  field  causes  a  net  force  on  the  other  elements  in  an 
array . 

The  farfield  pressure  distribution  can  be  computed 
from  Eq.  (A-3)  by  substituting  into  it  the  asymptotic  form  of 
hm(kr)  for  large  r,  which  is 

.  -m- 1  •< 

Lim  h  (kr)  =  — * e 

m  7  kr 

kr  -  ® 


Thus , 


Lim  p(r,\Ji  )  =  p  cv 


i(kr-u)t) 


F(cos  i!i)  , 


(A-6) 


1C  -+  oo 


where 


F  (cos  i!r)  = 


1  i  Tn[Pm_>^(cos  o  ^  “  Pm+1  (cos  a)]  Pm(cos  ili) 

^  - 2 - - - 

m=0  h' (ka)  +  -3— (kb)  ^h' (kb)[  j  '  (ka)n' (kb) -j  '  (kb)n' (ka)) 
mx  7  ipcv  7  m'  7LJmv  7  mv  7  Jmv  7  m  ^ 


This  expression  determines  the  farfield  characteristics  for  a 
single  piston. 

The  farfield  (or,  similarly,  the  nearfield)  sound  pres¬ 
sure,  p(r,0,§),  of  an  array  of  pistor.s  can  be  found  by  superposi¬ 
tion  of  the  fields  produced  by  the  individual  pistons.  The  piston 
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with  center  at  0j,  $j,  ^ j  on  the  transducer  has 
The  cosine  of  the  angle  between  9  ,  $  and  9j, 


velocity  Vje  *llt: . 
is  just 


cos  ijij  =  cos  9 


cos  0 j  +  sin  0  sin  9j  cos  ($  - 


v 


(A- 7) 


Superposition  of  the  pressure  fields  of  an  array  of  N 
pistons  gives 


i(kr  -  a)t) 

p(r,0,$)  =  pc  - - - -  D(9,§)  ,  (A-8) 


where 


N 

D(0  ,$)  =  y  Vj  F (cos  j) 
J=1 


The  farfield  intensity  distribution  I  is  p*  p,  where  the 
asterisk  denotes  the  complex  conjugate.  Therefore, 

1  =  fiF  \  I  D(e  >4)l  2  '  (A"9) 
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